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Abstract. This article presents a new method to calculate eigenvalues of right 
triangle billiards. Its efficiency is comparable to the boundary integral method and 
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insight into the statistical properties of the spectra. We analyse numerically the 
correlations in level sequences at high level numbers (> 10 5 ) for several examples 
of right triangle billiards. We find that the strength of the correlations is closely 
related to the genus of the invariant surface of the classical billiard flow. Surprisingly, 
the genus plays and important role on the quantum level also. Based on this 
observation a mechanism is discussed, which may explain the particular quantum- 
classical correspondence in right triangle billiards. Though this class of systems is 
rather small, it contains examples for integrable, pseudo integrable, and non integrable 
(ergodic, mixing) dynamics, so that the results might be relevant in a more general 
context. 
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1. Introduction 

Polygon billiards have been studied both classically and quantum mechanically for 
roughly twenty years now [1]. These systems are situated right on the borderline between 
integrability and chaos. They are usually divided into two classes: the rational polygon 
billiards where all vertex angles are rational multiples of n, and the irrational ones where 
at least one vertex angle is an irrational multiple of n. 

In the first case, there exist two constants of motion, so that one would expect 
integrability. However, due to singularities in the billiard flow, the invariant surface of 
the flow is not necessarily a torus (with genus g — 1), but may be of a more complicated 
topology (1 < g < oo). This produces a very complicated classical dynamics (see: 
[2, 3, 4] and references therein). The systems are called integrable if g = 1 and pseudo 
integrable [1] otherwise. 

In the second case (the irrational polygon billiards), there is no second constant 
of motion. These systems are typically ergodic [2] and probably weakly mixing [5, 6], 
though the Kolmogorov-Sinai entropy [7] is always zero. 

Quantum and semiclassical calculations have been performed from the very 
beginning [1, 8, 9, 10], but only recently [11, 12] it became possible to calculate 
sufficiently large level sequences at sufficiently high energies, such that correlation 
properties could be analysed directly. There are fundamental open questions: 

(i) Do the correlations in the spectra of polygon billiards eventually become stationary 
at sufficiently high energy? 

(ii) Are there families of polygon billiards with common statistical properties 
(universality)? 

(iii) What is the signature of classical pseudo integrability in the quantum spectrum 
(quantum-classical correspondence) ? 

On the one hand, there has been numerical evidence [12], that at very high energies 
the spectra of irrational triangle billiards are statistically similar to spectra taken from 
the Gaussian Orthogonal Ensemble (GOE). On the other hand, based on the numerical 
study of the spectra of several rational right triangle billiards, it was proposed that 
pseudo integrability implies so called "intermediate statistics" [11]. For the nearest 
neighbour distribution [13] this means: linear increase at small spacings (as in the GOE 
case) and exponential fall-off at large spacings (as for a random Poissonian sequence). 
Intermediate statistics has also been found in the context of disordered systems at the 
metal-insulator transition point [14, 15, 16], which might indicate some relationship be- 
tween both classes of systems. 

This paper is mainly concerned with question (iii). We consider the one-parameter 
family of right triangle billiards, labeled by the value of the smallest vertex angle 
< a < it/ 4. For this class, a secular equation is derived, which identifies the eigenvalues 
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as zeros of the determinant of a particular matrix K(E). Though the matrix is infinite, 
its elements are given explicitly by very simple expressions. This makes K(E) an ideal 
point of departure for numerical and analytical studies. 

The most obvious characteristic of rational polygon billiards is the genus g of the 
invariant surface of the classical Hamiltonian flow (the irrational polygon billiards can 
be included, setting g = oo). Hence we will investigate in detail the relation between 
g and the correlations in the quantum spectra. In the numerical part, level sequences 
are calculated at absolute level numbers > 10 5 for various examples of right triangle 
billiards. This provides valuable complementary information to recent results from Bo- 
gomolny et al. [11]. In the analytical part, the matrix K(E) itself is considered. Though 
K(E) is a pure quantum mechanical object, it is shown that g and 7 (which is closely 
related to g) play a crucial role for iterated mappings of the form \l/(n) = K n (E)ty(0). 
Based on this observation, a mechanism is proposed, which can explain the connection 
between the genus g and the correlation properties of the quantum spectrum. 

In section 2 a secular equation is derived for the calculation of the eigenvalues of 
right triangle billiards. It is used in section 3 to obtain and analyse the level spacing 
distributions for several right triangles. In section 4 we analyse the properties of the 
matrix K{E) itself, and we discuss the roles of the two classical parameters g and 7 in 
this context. The conclusions are presented in section 5. 

2. Secular equation 

Our point of departure is the observation, that any right triangle can be obtained from 
cutting an appropriate rectangle along its diagonal. This is used to derive a secular 
equation of drastically reduced dimension for the eigenvalues of the right triangle billiard. 

Let H Q be the Hamiltonian for the rectangle billiard with sides a and b. Fixing 
the length scale by: a 2 + b 2 = n 2 , the angle a : tana = b/a suffices to characterize the 
system completely. Choosing an arbitrary corner of the rectangle billiard as the origin 
of a Cartesian coordinate system, its eigenvalues and the corresponding eigenfunctions 
may be written as follows: 

1 / n 2 m 2 \ 

e( n ,m) =- — — + — ^— , n,m>\ (1) 

2 \cos 2 a; sin a / 

®nm(x, y) = -y= sm(^-nx^) sin(^my) . (2) 

Consider the total Hamiltonian H: 

H = H + V W , W = 6(^-^ , (3) 

where the potential ryW is used to cut the rectangle billiard into two congruent right 
triangle billiards (a similar cut potential, though in a different context, has been used 
in [17]). As i] increases from to 00, the spectrum of H changes from the spectrum 
of the rectangle billiard (1) to the doubly degenerated spectrum of the two triangle 
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billiards. For any 77, the Hamiltonian H is invariant under point reflection, so that 
the matrix representation of H in the eigenbasis of H is block diagonal. One block is 
spanned by the odd basis states {$ n m| n + m : odd} and the other by the even ones 
{$ nm | n + m : even}. Both blocks can be diagonalised independently, leading to the 
same sequence of eigenvalues, which causes the degeneracy mentioned above. 

In what follows we will work in the odd basis only. Let q = n + m and p = n — m, 
and order the states (2) with increasing q, and for equal q, with increasing p. Consider 
the subset of states with fixed q and p = — q + 2, . . . , q — 2 as one block. Then truncating 
the basis at a maximal q- value g max , one obtains M = (g max — l)/2 blocks with q — 1 
states in each block (note that q and p are odd). In total this gives N = (g^ ax — l)/4 
basis states. In this reordered basis, the matrix elements of W are given by: 

W qp . qV = / dx dy $ nm (x,y)$ n , m ,(x,y) 8 ( - - |) 

Jo Jo Vo b) (4) 

= \{K\V\ - \P'\) + S(q - \p'\) + 5(q' -\p\) + 5(q - q')} , 

where n = (q + p)/2, m = (q — p)/2 (and similarly for the primed indices). For 
given g max the truncated matrix has only two distinct eigenvalues: and M + 1, 

and the eigenspace of the latter has dimension M (in other words: rankfW^^] = M). 
All eigenvectors with eigenvalue M + 1 can be calculated explicitly, and after proper 
normalization we collect them (as column vectors) in the rectangular matrix V: 



Vk;qp 



: n q < k 



fc+i 

k ■ "'/ 

1 



VM+1 



k(k+l) 



n q = k 

n q > k , n p < k + 1 , (5) 



k 



, +1 : n q > k , n p — k -\- 1 
: n q > k , n p > k + 1 

where k = 1,...,M, n q = (q — l)/2 and n p = (\p\ + l)/2. The truncated total 
Hamiltonian may now be written as follows: 

H (N) = H (N) + ^ M + x) yyT _ (6) 

Dividing the Schrodinger equation (E — H^)^f = by E — Hq N \ one arrives after a few 
algebraic manipulations at the desired secular equation. It determines the eigenvalues 
of H W as the zeros of the following determinant: 

= det (l +r l K^ M \E)) , K {M \E) = (M + 1) V T . (7) 

E — H 

Taking the limit i] — > 00, the unit matrix in the first equation of (7) can be neglected, 
and one gets: 

det K {M \E) = 0. (8) 
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The advantage of this equation, is the reduced dimension M ^ N = MiM + 1). Such a 
reduction is typical for a boundary integral method (see for example [18]). The matrix 
elements of K^ M ^ are given by the following expression: 

KP = (M + 1)J2 £ f^g^ , ^:odd. (9) 

g=l p=- g +2 fc I 2 ' 2 J 

Being only interested in the zero eigenvalues of K^ M \E), any (symmetric) similarity 
transformation = L T K^L may be applied. The following choice for L simplifies 

the problem considerably: 

/ 1 -1 \ 



L = diag(l,...,l/M) 



-1 
1 



diag(l,..., v / M(M + l)) . 



(10) 

The resulting matrix K^ M \E) is defined in the same way as K(E) in the equations (11)- 
(13) below, but with the coefficients Dj = X^o^i- Only in the limit M — > oo, the 
expression for Dj simplifies to the formula (14), as can be shown using the partial frac- 
tion expansion of the cot function [19]. 

To summarize, the right triangle spectrum is calculated using the secular equation 
det[K(E)] = 0, where K(E) is constructed as follows: 

K{E)=K F {E) + K D {E) . (11) 

The matrix elements of K F are given by: 

[-^F]ij = — di j + \ — di—ij + di_ij + i , 

dij^d+j + d-;, d%= e _ g 2_ p 2 ±2qpcos2a > 

where the scaled energy e = E/{2 sin 2 2a) is used, and q — 2i + 1, and p = 2j — 1. Note, 
that q 2 +p 2 =F 2qp cos(2o;) = 2sin 2 (2o;) e[(q±p)/2, (q=fp)/2], thus q and p may still be 
regarded as auxiliary quantum numbers for the rectangle billiard H . The matrix K D 
is tridiagonal: 

[Knjjj = Dj + D j+1 , [K D \j, j+1 = [K D ] j+hj = -D 3 , (13) 

with the coefficients Dj, given by 

„ 7T sin7rcj / . , 

Dj = -, uj = sin 2a W 2e - p 2 . (14) 

2ui (cos huj + cos np cos 2a) 

Even though uj becomes imaginary for large values of p, the affected functions: sin and 
cos convert into sink and cosh, and finally the coefficient Dj remains real. Its asymptotic 
behaviour for large j is: Dj ~ 7r/(2|u;|). 
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This result is the basis for the analysis in section 3 and section 4. The infinite 
matrix K(E), as defined in (11)-(14), completely determines the spectrum of any right 
triangle billiard as the set of zeros of its determinant. For numerical purposes K(E) 
must be truncated (see below), but one may get important information also from an 
analysis of the infinite matrix K(E) itself. 

For numerical calculations (section 3), K(E) is truncated, keeping only those 
elements K^E) for which i, j < M. For meaningful results, M must be at least so large, 
that p^ ax > 2e, p max = 2M — 1 [see the definition of p above equation (13)]. Experience 
shows, that for accurate results (error less than 1% of the mean level spacing), one 
should increase the size of the matrix further by approximately 10%. The zeros of 
det[K(E)] are identified, calculating the smallest eigenvalue in magnitude as a function 
of E. Using a standard root bracketing algorithm [20] we find those points at which the 
smallest eigenvalue of K(E) passes the zero axis. The eigenvalues of K(E) are strictly 
decreasing functions of E, and this facilitates the root finding considerably. It allows to 
take rather large steps (of the order of the mean level distance), without running the 
risk to loose any roots. 

3. Level spacing distributions 

In the case of polygon billiards, the genus g of the invariant surface of the Hamiltonian 
flow is the most obvious parameter to characterize the classical dynamics [2]. Hence 
one may expect an influence of g on the level statistics of the corresponding quantum 
system. In this section we investigate numerically whether the level statistics show a 
systematic dependence on g. For several rational and one irrational right triangle, we 
calculate sequences of 10 4 levels starting at the absolute level number 10 5 (Weyl's law is 
used to determine the corresponding energy). Note that even in this energy region the 
level statistics are usually not stationary. This has been demonstrated in [12] for several 
examples of rational and irrational triangle billiards. This should be kept in mind in 
the discussion of the numerical results. 

In the case of right triangle billiards, there is another relevant parameter intimately 
related to g (see Appendix A). This is 7, the smallest integer such that 2a 7/^ e N (in 
the irrational case, we set 7 = 00). It is shown in Appendix A, that 7 is the smallest 
number of rhombuses which must be glued together to form the invariant surface of the 
billiard flow. Moreover we find, that g = int(7/2). Hence 7 implies a finer classification 
of the right triangle billiards than g does. 

Up to an irrelevant energy scale, all right triangles may be labeled and identified 
through their smallest vertex angle < a < 7r/4. For rational right triangles one may 
also use the pair of relatively prime integers p/q — a/ir. This is done in table 1, where 
all rational right triangles with g < 7 are arranged with increasing 7 in the vertical 
direction, and with increasing a in the horizontal direction. The parameter g a is taken 
from [11], where it is introduced as "arithmetical genus". The entries under-laid with 
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Table 1. All rational right triangle billiards with g < 7, referenced by their smallest 
vertex angle a/n — p/q, and ordered with respect to g and 7. The first two entries 
for g = g a = 1 are the only integrable cases. The shaded entries refer to those cases 
analysed in [11], and the gray-scale corresponds to the value of g a (introduced there) 
as indicated in the last column. 
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a gray shade have been analysed there. The gray-scale corresponds to different values 
of g a , as indicated in the last column. Further below, we will compare our results with 
those of [11]. 

We analyse the level statistics by means of the nearest neighbour spacing 
distribution P(s) (the spacings are normalized to unit mean). Recently the question 
has been raised, whether the rational right triangles show intermediate statistics (i.e. 
a linear increase at small s and exponential fall-off at large s). A simple analytical 
example is the so called "semi-Poisson" distribution [11]: 

P SP (s) = 4s e" 2s . (15) 

Here, P SP (s) is simply used as a conveniant reference to compare with. The following 
quantity is plotted in figure 1: 

A/ (s) = f ds' {P Num (V) - Psp(s')} • (16) 
Jo 

The theoretical curves in the figures 1(a) and (b) show the result for an infinite GOE 
spectrum, where P Num (s') is replaced by the corresponding level spacing distribution 
(taken from [21]). While figure 1(b) shows the raw numerical data curves for various 
right triangle billiards, figure 1 (a) shows the corresponding smoothed curves, in order to 
allow the identification of all the cases shown. For the smoothing, "natural smoothing 
splines" have been used, as provided in [22]. 

Let us first focus on the cases: p/q = 1/8,1/5,1/12,1/7,3/16 (which correspond 
to a successive increase of 7 from 4 to 8), and a/ir = (3 — v / 5)/4 (where 7 = 00). The 
AJ-curves for these cases are plotted in figure 1 (a) with a solid line and dashed lines of 
different dash lengths. Together with the results for the GOE and the Poisson ensemble 



Spectral properties of right triangles 



8 




Figure 1. Difference A J (s) of the integrated level spacing distribution to the scmi- 
Poisson case. AI (s) is plotted for various values of a as indicated in (a). The 
abbreviation "irr" refers to a = 7r(3 — \/5)/4. Whereas (b) shows the raw data, (a) 
shows the corresponding smoothed curves (details in the text). The GOE expectation 
(N — > oo limit) is plotted in (a) and (b) as a dotted line. 



(uncorrelated random sequence), they roughly span a one-parameter family of curves 
AI a (s). The parameter a may be called "correlation strength" and it may be calibrated, 
requiring that o = gives the Poisson result (its graph is plotted in figure 4), a — 1 
the GOE result, and o = 1/2 the semi-Poisson result. Note, that AI a (s) is introduced 
solely to facilitate the discussion of our results, so that it is not necessary to be more 
specific. 

The A/-curve for the irrational right triangle billiard comes closest to the GOE 
result. However, it remains almost in the middle between the semi-Poisson case and 
the GOE case (a > 3/4). Then follow the cases p/q = 3/16, 1/7, and 1/12, for which 
a decreases in approximately equal steps. The Ai-curve for p/q = 1/12 is closest to 
the semi-Poisson result (a ps 1/2). The last two curves with p/q = 1/5, and 1/8 tend 
slightly towards the Poisson result. They are so close to each other, that we would 
assign the same correlation strength to both of them (<r < 1/2). In all we may say, that 
the correlation strength a increases with increasing 7. 

Finally we included two more cases: p/q = 1/10 and 3/14. In figure 1(a) the 
respective AJ-curves are plotted with dotted lines. Thus we can compare the AJ-curves 
for the 1/5- and the 1/10-triangle (7 = 5), and the AJ-curves for the 1/7- and the 
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3/14-triangle (7 = 7). Both cases show, that even for right triangles with the same 
value for 7, the respective A/-curves may differ considerably. The relation between 7 
and the correlation strength is apparently not very strict (at least not in the energy 
range considered). 

In order to check, that our conclusions do not depend on the particular choice of the 
correlation measure, we repeated the analysis above, using the number variance S 2 (/) 
[13] instead of AI(s). The results were perfectly compatible, so that a more detailed 
discussion is omitted. 

In the numerical analysis presented here, we concentrate on right triangle billiards 
with small values for g and 7. The main reason is, that there is certainly an energy 
scale below which the quantum system cannot possibly "recognize", whether the two 
hypotenuse angles a and f3 are rational or not. Without any knowledge about this scale, 
small g triangles are probably the more reliable examples, for the study of quantum 
signatures of pseudo integrability. Hence the triangles for which we can compare our 
results with those of [11] are only a few: p/q = 1/8, 1/5, 1/12, and 1/7. 

The results obtained in [11] agree with those presented here, only up to a certain 
qualitative level. Beyond, we find that the correlation strength has decreased consid- 
erably in almost all cases. This may be due to the higher energy region considered 
here, which results in the rationality of the hypotenuse angles being more important. 
However, the Ai-curves of the first group of right triangles (with 4 < 7 < 6) changed 
much less then the others, such that the separation between both groups has decreased. 
It seems that this separation was decisive for the introduction of the arithmetical genus 
g a . The fact that this separation has become much smaller now, indicates that g a is 
probably not an appropriate alternative for g. 

According to the numerical results presented in this section, it is possible to order 
the right triangle billiards with respect to the strength of the correlations found in 
their spectra, which coincides with that of increasing 7. This finding confirms the 
general conjecture, that the genus of the invariant surface of the classical billiard flow 
determines the strength of the spectral correlations on the quantum level. Though the 
spectral correlations are apparently not stationary at currently accessible energies, the 
ordering seems to be energy independent, as long as the level sequences to be compared, 
start with the same absolute level number. 

4. The elliptic map 

In the first part of this section it is shown, that the parameters g and 7 (see Appendix 
A) associated with the classical dynamics of right triangle billiards, are important char- 
acteristics of the matrix K(E) itself. On the one hand this may be surprising, because 
K(E) arose from a pure quantum mechanical approach (see section 2), but on the other 
hand it is a strong indication for the importance of classical pseudo integrability on the 
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Figure 2. Portrait of the matrix K(E) for a = tt/5, and E = 1.26 x 10 4 . The 
gray-scale corresponds to the absolute value of the matrix elements. 



quantum level. In the second part, we present a tentative explanation for the depen- 
dence of the spectral statistics on g and 7. 

In figure 2 the matrix K(E) is portrayed for a typical case. The gray-scale 
corresponds to the absolute value of the matrix elements. It can be seen, that most 
of the matrix elements have vanishingly small absolute values. Large absolute values 
can be found only along the diagonal and the first off-diagonals, which are due to Kp, 
and on a "moon"-like structure due to Kp [see equations (12)-(14)]. The matrix elements 
become large, whenever the pair of integers is close to the zero-line of one 
of the two functions 



f± ( x ,y) = e — 4 (x 2 + y 2 ± 2xy cos 2a) , 



(17) 



where x,y are real, and positive, and e is the scaled energy as used in equation (12). 

The action of K{E) on a localized state may be described schematically by a double 
valued, symmetric map as shown in figure 3. The square in the middle represents the 
matrix K{E) (cf. figure 2). An initial state y localized at a given value a is mapped 
to yi = K(E)y Q localized at {ao,a_i,ai}, where a_i and a\ are the two solutions of 
the equation f±(a ,x) = 0, for x > 0. Hence, the map M associated with K(E) may 
be defined as follows: {a_i,ai} = Ma . Let us call it the "elliptic map". Due to 
f(x,y) = f(y,x), ao e Ma_i and ao G Ma\. Consequently, an orbit of such a map may 
be viewed as a doubly connected chain: 



a- 3 - 



-a-2- 



-a_i- 



-11- 



-02- 



-a 3 



(18) 



According to that picture, the n-fold image y n = K(E) n y has localization peaks 




Figure 3. Schematic representation of the mapping y\ = K(E)y (details see text). 



at the positions {a_ n , . . . ,a , . . . , a n }. Surprisingly, M is isomorphic to the following, 
extremely simple map: 

Wti = V?n ± 2a , (19) 

where the result should be taken modulo ir, such that it remains in the interval 
[ — 7r/2,7r/2). This can be seen, using the following parametrisation of the curve 
f±(x,y) = 0: 

pUv^f Sgn(f ~ 2a + f\ COS(v? - 2a) V <pe[-«/2,«/2) .(20) 

\yj V c os(v^) / 

Replacing y by an arbitrary point a n of the map M, one gets the corresponding pair 
conjugated angles: cos(±y? n ) = a n . Replacing x by a n one finds that ±</? n must be 
mapped to ±</? n — 2a (mod n). It follows, that 

a n = cos(±y? n ) 

Ma n = {a n _i, a ra+ i} = {cos[±(y? n - 2a)], cos[±(</?„ + 2a)]} . (21) 

From equation (19) it follows that any orbit is restricted to a set of 7 points, where 7 
is the smallest integer such that 2a 7/71 e N. It is the same 7, which is introduced in 
Appendix A as the number of rhombuses forming the invariant surface of the billiard 
flow. Furthermore, the periodicity of the map M is hit (7/2) which is just the genus of 
that invariant surface. 



Here in the second part of this section we discuss a mechanism which can explain 
the correspondence between the correlation properties of the quantum spectrum and 
the classical parameter 7. The line of reasoning is as follows: 

1 The correlation properties of the triangle spectrum at a given energy E are closely 
related to the correlation properties of the eigenvalues of K(E) in the vicinity of 
zero. 

2 At sufficiently high energy, K D (11) can be considered as a random tridiagonal 
matrix with eigenstates which are typically localized. 




s 



Figure 4. Difference of the integrated level spacing distribution to the semi-Poisson 
case. In (a) AI (s) is plotted for the rational right triangle billiard with a/ir = 1/8, in 
(b) for the irrational one with a/ir = (3 — \/5)/4. The solid lines show the result for 
neighbored eigenvalues of the matrix K(E), the dotted lines show the result for the 
triangle spectrum. The dashed lines show the theoretical curves for the Poisson case 
(long dashed lines) and for the GOE case (short dashed lines). 



3 The matrix K F (11) has such a form, that repeated multiplications of an initially 
localized state with K(E), produce an increasing number of copies at different 
positions. The positions are given by the elliptic map M. 

4 If a is rational, all orbits of the elliptic map are periodic with period g and restricted 
to 7 points. This leads to an approximate foliation of the Hilbert space into weakly 
coupled subspaces. For any irrational a, the elliptic map is ergodic, and all basis 
states of the matrix K(E) are strongly coupled. 

Point 3 has been treated in the first part of this section. The remaining statements are 
discussed below. 



Correlation properties of the eigenvalues of the matrix K{E) 

According to the secular equation derived in section 2, the triangle eigenvalues are given 
by those energies, at which at least one eigenvalue of K(E) becomes zero. Therefore 
it seems plausible, that the correlation properties of the eigenvalues of K{E) close to 
zero and the triangle eigenvalues are related. To verify this we calculate the spacing 



Spectral properties of right triangles 13 




50 100 150 200 250 -2 -1 1 2 

X10- 4 



Figure 5. The eigenstates of Kd(E) with eigenvalues close to zero for a typical case. 
On the left, the absolute value of the eigenvector coefficients, plotted as a function of 
i, the index for the basis of K{E). On the right, the corresponding eigenvalues plotted 
in a bar graph, a = n/5, E = 1.5 x 10 5 . 



distribution for those two neighbouring eigenvalues which have opposite signs (without 
unfolding). The distribution is obtained from 10 4 spacings, taken at equi-distant 
energies, with the step size adjusted to the mean level spacing of the corresponding 
triangle spectrum. 

In figure 4 we compare the results. In the rational case (a) as well as in the irrational 
case (b), the A/(s)-curves for the eigenvalue pairs of K(E) and for the triangle spectra 
differ remarkably, though in (b) the agreement is somewhat better. However, at least 
qualitatively, the results are as expected: In the rational case, figure 4(a), the eigenvalue 
statistics for K{E) show indeed very weak correlations, much weaker even than the 
corresponding triangle spectrum. This can be seen from the A/(s)-curve which clearly 
tends towards the Poisson result (note also the behaviour at large s). In the irrational 
case, figure 4(b), both curves show relatively strong correlations. 

Localization of the eigenstates of Kr> 

The matrix Ke> is constructed in a simple manner from the coefficients Dj [see 
equations (13) and (14)]. Those oscillate as functions of j, the index of the basis of 
K(E), more and more rapidly while p 2 approaches 2e. From a statistical point of view, 
it then seems permissible to replace the arguments of the functions: sin and cos by 
appropriate random variables. Even though the statistical properties of the matrix 
elements are very complicated, one may expect Anderson localization [23]. 

In figure 5 we show for a typical case, a series of eigenstates of K D (E) ordered 
by their respective eigenvalues. Only those states with eigenvalues in a small interval 
around zero are shown. Many eigenstates are apparently localized. However, others are 
not, and spread over a wide range of basis states. Usually those fluctuate only weakly 
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and slowly and their eigenvalues decrease very slowly with energy (not shown). Their 
role is still unclear, and will be the subject of future studies. 

Foliation of the Hilbert space 

Acting repeatedly with K{E) on an initially localized state y , the first g images will 
spread and localize at points {ao, . . . , a 7 _i} (here a 7 _i is identical to a_i), as discussed in 
the first part of this section. Then, due to the periodicity of the elliptic map, subsequent 
images spread only slowly away from these points. In the ideal case the spreading would 
stop due to Anderson localization, giving rise to an invariant subspace. In the same 
way, an initial state localized in a different part of the Hilbert space, would lead to 
another invariant subspace, and so on - until possibly the whole Hilbert space would 
have been decomposed into invariant subspaces. In the real system, such a foliation of 
the Hilbert space occurs only approximately, and the subspaces become weakly coupled. 
Nevertheless one may expect, that correlations are to some extent suppressed due to 
this mechanism. 

For irrational a, where the elliptic map is ergodic, an initially localized state will 
spread out [by repeated multiplication with K(E)] into the whole Hilbert space. No 
foliation of the Hilbert space can occur, and one should expect correlations of similar 
strength as in the GOE case. 

5. Conclusions 

We derived a new kind of secular equation for the determination of the spectra of 
right triangle billiards. It involves the diagonalisation of the matrix K(E) which has 
a particularly simple and transparent structure. Based on this equation we calculated 
spectra at level numbers > 10 5 for various examples of right triangle billiards, which 
shows the efficiency of the new method. 

We found a clear correspondence between the genus g (or the related parameter 7) 
of the invariant surface of the classical billiard flow and the strength of the correlations in 
the quantum spectrum. While for small g the spectral statistics is close to semi-Poisson 
(with a slight tendency towards Poisson), it approaches the GOE statistics when g is 
increased. Our numerical results together with similar studies [11, 12] suggest that the 
spectral correlations are not stationary at currently accessible energies, but that the 
ordering with increasing correlation strength and its correspondence to g is conserved. 

In the second part of the paper, we found that the classical parameters g and 7 are 
characteristic quantities for the matrix K{E) itself. For rational right triangle billiards, 
where g is finite, one gets an approximate foliation of the Hilbert space into invariant 
subspaces. The size of the subspaces scales with 7. Based on this observation, we dis- 
cussed a mechanism which can explain the influence of g and 7 on the level statistics of 
right triangle billiards. 
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The definition of 7 in Appendix A can be generalized to arbitrary polygon billiards 
as follows: 7 is the smallest number of identical polygons which must be glued together 
to form the invariant surface of the billiard flow. In this general case, g and 7 must 
possibly be considered as independent parameters. Which of them is then more relevant 
for the correlations in the quantum spectrum? Considering the roles of g and 7 in the 
description of the matrix K(E), it seems that this is 7 (the size of the approximate 
invariant subspaces of K{E) scales with 7). 
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Appendix A. The invariant surface for the classical billiard flow 

There is an elegant way to represent a trajectory moving in a polygon billiard, which 
is particularly useful to construct the invariant surface of the classical billiard flow. It 
consists in drawing the trajectory as a straight line, and reflecting the billiard (instead 
of the trajectory) each time the boundary is hit [2]. In the case of rational polygons all 
possible trajectories can produce only a finite number of differently oriented copies of 
the original polygon. Then there is a general recipe of how to glue these copies together, 
in order to obtain the invariant surface. 

For rational right triangles one may follow a more explicit construction scheme, 
which leads to a particularly simple invariant surface; the "rosette" . It is constructed as 
follows: Start with a right triangle as depicted in figure Al(a). Reflecting the triangle on 
the side AC, the image on the side BC and that image again on the side AC, produces 
a rhombus [see figure Al(b) and (c)]. This rhombus is rotated around point A by the 
angle 2a (counter-clockwise) at each step. Stop one step before arriving at the original 
rhombus or its point reflected image [see figure Al(b)]. Note, that the resulting surface 
may wind several times around A. As it is shown below, the resulting surface can be 
closed by identifying open edges with one another, which gives the invariant surface. 

The number of rotations by 2a is just the smallest integer 7 such that 2a r y/ir G N. 
We could equally well rotate the rhombus step-wise by 2(3 around point B [see 
figure Al(b) and (c)], resulting in a different representation of the invariant surface. 
However, as shown below, the number of rotations (or rhombuses) necessary to close 
the invariant surface is the same. 
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Figure Al. (a) Typical right triangle, (b) Invariant surface of the billiard flow for 
the 1/8-triangle. The dotted rhombus does not form part of the invariant surface, (c) 
Invariant surface of the billiard flow for the 1/7-triangle. In (b) and (c): the initial 
rhombuses, with which the rosette constructions are started are filled. Vertices are 
labeled by capital letters, edges by lower case letters. However, only such edges or 
vertices are labeled, which must be identified with one another, to obtain the invariant 
surface. 



Consider the first line of (A.l). As 2p < q it follows 7 — Z > 0. But 7' is the smallest 
integer such that 27' p'/q' G N, hence: 7' < 7. The same argument applied to the 
second line of (A.l) shows: 7 < 7'. Therefore: 7 = 7'. 

It remains, to prove that the procedure above gives indeed a representation of 
the invariant surface, and to calculate its genus g. To this end we show that all free 
edges of the rosette can be identified with one another. Then the rhombuses define 
a triangularisation of the invariant surface, and counting all faces F, edges E and 
vertices V of the triangularisation, we obtain the genus via the Euler characteristic [2]: 



The edges may be divided into inner edges, which are connected to the center of 
the rosette, and outer edges, which are not connected to the center. Let us label both 
groups counter-clockwise by ei,. . . and e' lv . . respectively, beginning with the lower edges 
of the initial rhombus (see figure Al, but note that the labels shown there are different, 
and used only to identify different edges in the representation of the invariant surface). 

Let us first discuss the case, where q is odd. Then 7 = q and the rosette winds p 



Proof: Let p/q = a/ir and p'/q' — fl/ir with p, q and p', q' relatively prime. 




(A.l) 



g = l- X /2,x = V-E + F. 
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times around its center before the last inner edge can be identified with the first one [see 
figure Al(c)]. In order to identify all outer edges ei,. . . ,e 27 pairwise with one another, 
observe that for j :odd, a trajectory leaving the surface crossing e J+ 3 would enter a 
triangle which is the parallel translation of the triangle with the hypotenuse e^. By 
consequence, both edges can be identified. Hence one may identify the following edges: 
ei = e 4 , e 3 = e 6 , . . . , e 27 _ 3 = e 27 , and there are only two open outer edges left: e 2 and 
e 2 7 -i which can be identified with one another on the same grounds. 

The vertices in the representation of the invariant surface have to be identified 
taking into account that edges identified previously have the same initial and end points 
(the triangle connected to the edge defines an orientation). In this manner it is shown 
that for q :odd, the rosette (invariant surface) has 7 faces, 27 edges and 3 vertices. 
Hence X = 3 — 7 and g — (7 — l)/2. 

If q is even, then p (relatively prime) is odd and the rosette winds p/2 times around 
its center [see figure Al(b)] which means, that we have also two open inner edges: e[ 
and e 7+1 . Identifying outer edges as explained above, leaves two outer edges open: e 2 
and e 27 _i. In this case we identify e' : with e 27 _i and e 2 with e' +1 , which again closes 
the invariant surface. Note that due to the identification of outer edges with inner 
edges, the central point of the rosette must be identified with the outermost points. The 
remaining points, must be identified as one single vertex B, if 7 is even [see figure Al(b)], 
otherwise they constitute two vertices B and B' (this case is not shown). Hence the 
rosette (invariant surface) has 7 faces, 27 edges, and 2 vertices if 7 is even and 3 vertices 
if 7 is odd. This gives x = 2 — 7 =^ g = 7/2 and x = 3 — 7 =^ g = (7 — l)/2 respectively. 

In all: g — (7 — l)/2 if 7 :odd, and g = 7/2 if 7 :even. Hence g = int(7/2). 
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